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Abstract 

We study the systematics of neutron star cooling curves with three represen- 
tative masses from the most populated interval of the estimated mass distribu- 
tion for compact objects. The cooling simulations are made in the framework 
of the nuclear medium cooling (NMC) scenario using different combinations of 
possible nucleon-nucleon pairing gaps. Possible heating or enhanced cooling 
mechanisms in the crust are not considered. We define a constraint on the 
highest possible temperatures for a given age of young neutron stars and show 
that this limits the freedom of modeling pairing gaps and crust properties. 
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Studies of neutron star cooling evolution become very actual due to the presently known 
surface temperature and age data provided by X-ray observatories such as CHANDRA, 
XMM Newton and from the ROSAT catalogue [1]. These new data open a wide perspective 
for nuclear astrophysics for which up to now the knowledge of internal structure of the 
compact stars and the properties of stellar matter under extreme conditions remain central 
problems. Theoretical models and hypotheses about the equation of state of high density 
matter provide different alternatives for the explanation of the same set of observational 
temperature - age (TA) data points, when additional constraints are not provided. In this 
work we point out an existing correlation between the crust model and cooling behaviour of 
light neutron stars, which has a selective power in combination of TA data with the mass 
spectrum of neutron stars. 

In our recent investigations of the cooling evolution of neutron stars (NS) we have adopted 
the so called nuclear medium cooling (NMC) scenario [2], which goes beyond the minimal 
cooling scenario [3], where in-medium modifications of cooling regulators by definition have 
been disregarded. Both approaches agree in the philosophy that such very effective cool- 
ing mechanism like the direct Urea process should not occur in typical NS. In these ap- 
proaches the main cooling process is the modified Urea process, which in our NMC scenario 
also includes the in-medium softening of the pion propagator [4]. Earlier investigations 
within this cooling scenario [5,6] have chosen the crust model as a simplified Tsuruta law 
T s Tsur = (10 T in ) 2/3 . Although it is shown in Ref. [7] that the cooling evolution could be 
essentially affected by the inclusion of internal heating, nevertheless the latter is expected to 
be important for late time evolution and will not affect the results of this paper. It has been 
omitted in the present cooling scenario. For more recent reviews on the cooling scenarios 
see [8-10]. 

The cooling simulations presented in this work are based on a code with a number of 
improved inputs concerning the heat conductivity, the nucleon-nucleon pairing gaps and a 
new model of the neutron star crust and envelope. These models are basically taken from 
the recent calculations of Ref. [11], where the amount of light elements in the crust and the 
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influence of the magnetic field have been taken into account. 

The choice of the structure of crust and envelope of the compact object becomes a central 
question due to the direct connection of the surface temperature with the observations. 

Since simultaneous measurements of the surface temperature and either the mass or 
the radius of neutron stars are absent, one needs to relate the known TA data with other 
observational information in order to discriminate between different models. 

One quantitative approach has been suggested in [12,13] and developed recently for the 
NMC scenario [14]. It uses the statistical method of population syntheses in combination 
with models for the cooling evolution and the mass distribution of isolated objects to derive 
a Log N - Log S dependence to be compared with observational data. This allows for 
an additional selection among those scenarios which have successfully passed the TA test. 
Nevertheless, in the routine of the Log N - Log S test the early stages of the time evolution 
are not included, but only the tracks of relatively old objects (t <; 3 • 10 4 yr). This period is 
already in the transition from the neutrino to the photon cooling era. We will show in this 
work that the cooling behavior of the young objects is strongly affected by crust properties, 
which however do not alter the Log N -Log S distribution. 

In the spirit of the Log N - Log S test we suggest an additional condition to the TA test, 
which requires young neutron stars (10 3 < t/yr < 3 • 10 4 ) with low masses ~ 1.1M not 
to have temperatures (i.e. brightness = luminosity / 4 n distance 2 , where absorption could 
be neglected for hot sources since they have a hard spectrum) exceeding those of presently 
known objects of the same age. We will therefore call it a brightness constraint (BC). 

Therefore, it is unlikely that objects with a given age are hotter than those already 
observed. As a consequence, cooling models which would predict such objects should be 
rejected. 

Our approach is based on the following assumptions made on the basis of current obser- 
vational data: 

• All objects with temperatures exceeding T— 4 • 10 6 K are potentially observable. 
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The probability to find an object increases with its luminosity due to the increase 
of the visibility volume and reduction of the absorption of X-rays by the interstellar 
medium (see [15]). The absorption falls by several orders of magnitude with increase 
of the photon energy from 0.1 to ~ 10 keV. One can estimate the lower limit of the 
temperature of an object at the distance of 10 kpc or even more distant and behind 
clouds of the interstellar medium with maximal absorption, say for hydrogen column 
densities nn ~ 2 • 10 22 cm -2 (corresponding to galactic coordinates 1=30, b=0), which 
would be observed with instruments like ROSAT, with a lower limit for the brightness 
given by a few 0.01 cts/s. Using such an estimate NASA's HEASARC tools web- 
page 1 , one obtains a limiting temperature T = 4 • 10 6 K, close to the temperature of 
Crab pulsar (2.2 • 10 6 K). Therefore, in a good approximation the temperature of the 
observed object with highest temperature could be considered as a border of BC. 

To specify the BC as a test we have made further assumptions: 

• The distribution of young objects is a step-like function of temperature for higher tem- 
peratures. 

• The uncertainty of the BC border is of the same order as the error bars of the measured 
temperatures at the corresponding age. 

We assume that with sufficient observational data it is possible to measure the temper- 
ature distribution of the young objects. Even if the population of young objects is small 
(N ~ 100 [16]) one can not conclude that the probability to find an object with arbitrary 
high temperature is finite. The distribution is likely a step like or very stiffly falling func- 
tion of temperature for higher temperatures, because the star can not have arbitrarily high 
temperature from the beginning of the evolution and the mass distribution of the compact 
objects are also a step like function for the small masses (see [14]), which are expected to 



http:/ /heasarc. gsfc.nasa.gov / docs / tools.html 
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populate the domain of slower coolers. 

The statistical analysis to define the BC border is complicated by the small total number 
of potentially observable objects. It is assumed to be ~ 0.1*N [17], because all known objects 
are colder than 4 • 10 6 K [18]. 

Therefore we leave a corridor between the upper limits of the observed data points of 
higher temperatures at a given age and the BC border on the 2a level of the error bar 
of measured temperature where it was possible and put 5 log 10 (T/K) ~ (1/N)* ( the 
whole expected temperature interval of young objects)/(the probability to be observed) 
~ 0.4 * 0.01/0.1 = 0.04 (for Crab), which is of the same order as the error bars of measured 
temperatures. 

This last estimation is based on the assumption that the probability for the existence of 
an object is equally distributed for all available temperatures. 

In order to apply this simple BC test to already published [2] scenarios we choose the 
strategy to follow the cooling evolution of objects with the representative values M = 
1.1, 1.21, and 1.41M corresponding to the most populated bins in the mass spectrum of 
Fig. 1 in [14]. The basic idea behind this mass spectrum of NS is to use HIPPARCOS data 
on massive stars around the sun as the mass distribution of progenitors in conjunction with 
the calculations by [19] for a population synthesis of nearby NS. 

The choice of the above representative masses for the BC test is justified by the results 
of [2] and [3] where it was demonstrated that all objects outside the chosen interval of the 
masses 1.1 -j- 1.4M are cooler than those inside this interval. Even if this systematics is 
not too strong, this statement is based on the set of cooling simulations [2] and is applicable 
to young objects and can be rephrased as follows: For a given age the brightest objects are 
from the most populated mass interval. 

The description of the cooling evolution is mainly given by the equation combining the 
energy balance and thermal energy transport [20] 




(1) 
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where and A = — |ln(l — 2m/r) are metric coefficients. The heat conductivity k, the 
total neutrino emissivity e u and the total specific heat cy/n are given as the sum of the 
corresponding partial contributions defined for density profiles n(r) of the constituents of 
the matter under the conditions of the actual temperature profile T(r,t). The mass of 
the star is the accumulated mass below the surface, M = m(r = R), which together with 
the gravitational potential 0(r) can be determined by Oppenheimer-Volkov equations (see 
[6,20]), where the energy density profile e = e(r) and the pressure profile p = p(r) are 
defined by the condition of hydrodynamical equilibrium. The boundary condition for the 
solution of (1) reads T(r = Ri n ,t) = T in (t). The energy flux from the surface has also a 
contribution from photons L 7 = 4iiaR 2 T^, which is governing the cooling, when the inner 
crust temperature T in falls down to 10 8 K. 

For different crust models we used in our recent calculations the T s — T in relations [2,21,22] 
shown in Fig. I. The parameter rj ~ AM L /M is a measure for the thickness of the light 
element layer [23], which is related to the pressure at the bottom of the light element 
envelope. Thus these borders of the acceptable T s — T in relations can be denoted as heavy 
element rj = 4 • 1CT 16 (further denoted as crust model (E) corresponding to the notation in 
Ref. [3] ) and light-element rj — 4 • 10" 8 crust models, respectively. 

From the discussion we omitted the possible heating in the crust to make our focus on 
the general aspect of the choice of crust. In principle the heating could also be included 
without qualitative changes in the argumentation we use here, because it will change the 
cooling evolution only in the photon era, but not the evolution of young objects [24,25]. 

The surface temperature of NS with thick light element crust is higher during the neutrino 
cooling era and it shows slower cooling for young objects in contrast to the case with a heavy 
element rich crust (E) as shown in works [2] and [3]. 

For the cooling simulation done in Refs. [2,21] we follow the idea of light element decay, 
which is in coherence with Ref. [3] . The assumption is that during the time evolution of the 
envelopes chemical composition the mass fraction of the envelope consisting of light elements 
decays. The time dependence can be described by an exponential 
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AM L (t) = e^ T AM L (0) (2) 

where AM L (0) is the initial mass of light elements and r is the decay rate. This decay 
could be due to the pulsar mechanism which injects light elements into the magnetosphere 
or due to nuclear reactions which convert these elements into heavy ones [26,27]. This rather 
complicated picture we had modeled with a fit between the light and heavy element crust 
models corresponding to a slow decay of light elements (model (C)), whereas the heavy 
element crust model is denoted as fast decay (model (E)), see Fig. 1. 

For the demonstration of the results given in Fig. 2 we reconsider the same cooling 
scenarios as in [14] using the notation given in Table 1 of that paper. We exclude only 
models V and II from our present discussion, since they are not in our aimed class of 
models. The remaining models selected are defined in Table I by the choice of the nucleon 
pairing gaps: (A) models I, IX - gaps from [28]; (B) models III, VI, IV, VII, VIII - gaps 
from [11]. For all models the 3 P 2 neutron gap is suppressed by 0.1; see Refs. [2], [14] and 
[21] for more details. 

Each of the six models is calculated with both crust models (C) and (E). The results 
are displayed in the six panels of Fig. 2, grouped into two columns. The left one shows the 
models which successfully passed the TA test, but only for crust model (E) fulfill the BC 
for younger objects to not have higher temperatures than those given by the observations. 
In the right column all models pass the BC test, while those with crust (E) failed the TA 
test leaving some observed points out of explanation. Note that in Fig. 2 we have not 
shown cooling curves for neutron stars with masses exceeding 1.41 M which correspond to 
intermediate or fast coolers since we are focused here on the discussion of young and slow 
coolers. We want to stress, however, that all models in Fig. 2 can pass the TA test if the 
crust model is chosen suitably. 

In the upper two panels both models (IV,VI) and IX are calculated without possible 
7r-condensation in contrast to the models III and I, where the condensation is possible. 
Therefore, the difference in these corresponding plots is only in the behavior of the configu- 
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ration with M — 1.41 M for which 7r-condensation occurs. The difference between the four 
upper panels (III, IV) and (I, IX) comes from the difference in the nucleon-nucleon pairing 
gaps. Only the latter two models are microscopically justified [28]. However, without addi- 
tional suppression of the neutron 3P2 gap, necessary to reduce the enhanced cooling by the 
neutron pair breaking and pair formation processes, the models would be in disagreement 
with the TA test (see [21]). The physical reason of this suppression could be the medium- 
modification of the spin-orbit interaction in neutron matter [29] . Thus the idea of slow decay 
of light elements in crust model (C) could be considered as a more consistent suggestion of 
the crust-envelope model. 

Summarizing our discussion we can conclude the following. The application of the NMC 
scenario for the simulation of the cooling evolution of neutron stars in comparison with the 
existing observational data shows that 

- the TA test of cooling scenarios in conjunction with the BC test can be selective for 
the discrimination between the crust models, when nuceon pairing is already chosen, 
or vice versa; 

- the TA test of the cooling scenario even with the improvement of the BC test is not 
sufficient to make a final conclusion about the validity of neither the crust model nor 
the nuclear superfluidity; 

- since the Log N - Log S test essentially depends on the cooling behaviour of objects 
older than 3 • 10 4 yr it does not interfere with the results of the BC test. 

All three tests are necessary but not sufficient for the final selection. With a special choice 
of the crust model it is possible to change the results of the TA test while the LogN -LogS 
test will remain unaffected. As it is discussed in the Ref. [30] there is another accelerator of 
the cooling of low-mass neutron stars, the direct Urea process of neutrino emission allowed 
in the mantle of a neutron star near the crust-core interface, due to inhomogeneous nuclear 
structures. Such possibilities will make the crust model of slow decay of light elements more 
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preferable. 

In conclusion, the present work is a contribution to the development of general testing 
schemes for models of compact star cooling evolution [31] and strongly interacting matter 
at high-densities [32] using constraints from compact star observations. 
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TABLES 



Class of model 


Models 
with 7r-condensate 


Models 
without 7r-condensate 


Gaps 


A 


I 


IX 


[28]; n- 3 P 2 *0.1 


B 


III 


IV & VI 2 


[11]; n- 3 P 2 *0.1 


B' 




VII 


[11]; p- 1 S *0.5; n- 3 P 2 *0.1 


B" 


VIII 




[11]; p- 1 S *0.2; n- 1 S *0.5; n- 3 P 2 *0.1 



TABLE I. Classification of the models by the possible transition to it- condensate and the 
choice of proton (p- 1 So) and neutron (n- 1 So, n- 3 P 2 ) pairing gaps with corresponding suppression 
factors. 



2 The model VI is the same as model IV, which was already calculated with crust (E) in Ref. [2] 
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FIGURES 
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FIG. 1. (Color online) The relation between the inner crust temperature and the surface tem- 
perature for different models. Dash-dotted curves indicate boundaries of the uncertainty band. 
Notations of lines are determined in the legend. For more details see [2,3] and [23]. 
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FIG. 2. (Color online) Cooling evolution for representative NS from the most populated mass 
bins M = 1.1, 1.21, IAIMq according to the NMC scenario. The classification of models is taken 
from Ref. [14]. The cooling curves for slow decay (C) are shown with thick lines and those for 
fast decay (E) with thin lines. The data points correspond to Fig. 1 of Ref. [3]. According to the 
brightness constraint applied for young objects the shaded regions should not be populated and 
cooling scenarios entering there shall be rejected. 



15 



